Power-Law NLED-Based Magnetic Universe Can Mimic Phantom Behavior 
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We study the cosmic dynamics of a magnetic universe supported by non-linear electrodynamics 
(NLED) Lagrangeans that are proportional to powers of the electromagnetic invariant oc i^^/f^""') 
(m is an overall constant). For simplicity we focus in the case when F depends on the magnetic 
field alone, a case dubbed in the bibliography as "magnetic universe". Our results demonstrate 
that, depending on the values of the free parameter m, the magnetic field can mimic phantom field 
behavior, an effect previously found in other contexts. It is demonstrated that, since there are 
found equilibrium points in the phase space of these models that can be associated with magnetic- 
dominated past and future attractors, a combination of positive and negative powers of F may lead 
to interesting cosmological behavior. In particular, a cosmological scenario where the universe might 
evolve from a past NLED-driven (non-inflationary) state into a future (late-time) - also NLED-driven 
- inflationary stage, transiting through a matter-dominated solution, is envisioned. The impact of 
braneworld gravity on the dynamics driven by such NLED Lagrangeans is also investigated. It is 
demonstrated that, due to phantom property at late times, the non-linear electromagnetic effects 
may play an important role in deciding the fate of the cosmic evolution. Randall- Sundrum brane 
effects, in particular, modify the nature of the starting point of the cosmic evolution, as well as the 
fate of the expansion: both, the big bang singularity and the big rip event - inherent in general 
relativity with phantom matter source -, may be avoided. 

PACS numbers: 04.20.-q, 98.80.-k, 98.62.En, 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 



Studying the equations of the non-linear electrodynam- 
ics (NLED) is an attractive subject of research in general 
relativity (GR) thanks to the fact that such quantum 
phenomena as vacuum polarization can be implemented 
in a classical model through their impact on the prop- 
erties of the background space-time. Exact solutions of 
the Einstein's field equations coupled to NLED may hint 
at the relevance of the non-linear effects in strong grav- 
itational and magnetic fields. It has been speculated, in 
particular, that very strong electromagnetic fields might 
help avoiding the occurrence of space-time singularities 
in the cosmological context [l[ . The impact of very strong 
electromagnetic fields (and of the NLED effects) regard- 
ing the causality issue in cosmology is also of relevance 
A different subject of research within the cosmolog- 
ical setting in GR coupled to NLED, is related with the 
chance for the NLED field to fuel primordial inflation. 
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The cosmological inflationary scenario was proposed for 
the first time in the Ref. Q where, in order to over- 
come several problems of the standard model - such as 
the flatness and the horizon problems, amongst others -, 
it was anticipated that a self-interacting scalar field with 
a particular form of the self-interaction potential as the 
source of the Einstein's field equations, might cause the 
universe to expand in an inflationary (super-accelerated) 
fashion. The inflationary paradigm has been supported 
by the observational evidence In the reference 
for instance, an anisotropic cosmological model coupled 
to Born-Infeld NLED was explored. It was found that 
this model might explain early time infiation . Thanks 
to the infiationary stage, the initially anisotropic universe 
might eventually isotropize. An alternative non-isotropic 
model where, additionally to the Born-Infeld NLED field, 
a cosmological constant was added to the Einstien's field 
equations, has been studied in Q- Yang-Mills cosmol- 
ogy with a non-Abelian Born-Infeld action has been also 
considered 

Magnetic universes (vanishing electric component) 
have been investigated within the context of NLED the- 
ories given by the Lagrangean density Li — —jF + 
aF'^ + PG^ (see below) @. The non- linear term oc F^ 
may cause the universe to bounce thus avoiding the ini- 
tial (big-bang) singularity. However, in Ref.^lH] it was 
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demonstrated that the inclusion of a non-vanishing elec- 
tric component E where ~ B^, removes the bounce 
and, in consequence, the universe starts its evolution in 
a singular state. Lagrangeans with inverse powers of the 
electromagnetic field F are interesting because the non- 
linear electromagnetic effects might become important 
not only at early times in the cosmic evolution, but also 
at late times. Actually, models with Lagrangean density 
of the form L2 = —■jF— ^ may account for the late-time 
stage of accelerated expansion of magnetic universes [llj . 
In the later case, if a non- vanishing electric component E 
is considered, accelerated expansion is allowed only when 
E^ < 3^2 [g. 

Even if the above - very simplified - models of non- 
linear electrodynamics coupled to general relativistic cos- 
mology, describe hypothetical systems reminiscent of the 
fields in the real world, these models comprise interest- 
ing dynamical behavior that is worthy of independent 
investigation. Through the study of the equivalent phase 
spaces, for instance, it may be revealed how generic can 
be the occurrence of a bounce instead of the big-bang 
singularity, or, also, how the dynamics of these models is 
modified by the brane effects, etc. Aim of the present in- 
vestigation is to explore the asymptotic properties of cos- 
mological models where general relativity is coupled to 
NLED given by Lagrangeans that contain powers of the 
electromagnetic invariant F (Li and L2 above are par- 
ticular cases). To simplify the study we shall consider 
magnetic universes (vanishing electric component). In 
spite of their simplicity these models contain interesting 
physics. We will rely on the use of the dynamical systems 
tools with the hope that such relevant concepts as past 
and/or future attractors, or saddle equilibrium points, 
could be correlated with generic cosmological behavior. 
Our goal will be to write the cosmological (Einstein's) 
equations in the form of an autonomous system of or- 
dinary differential equations. Although these equations 
are non-linear, we will expand them in the neighborhood 
of the equilibrium (critical) points up to the linear ap- 
proximation. Evaluating the sign of the real parts of the 
eigenvalues of the corresponding linearization matrices 
will allow us to judge about the stability of the critical 
points, which are correlated with asymptotic properties 
of the original system of equations. 

The paper has been organized in the following form. In 
section II we expose the basic mathematical and physical 
aspects of NLED-based cosmological models. The cos- 
mic dynamics of models based on Lagrangeans that are 
power-law in the electromagnetic invariant F, is investi- 
gated in section III. For simplicity we concentrate in the 
so called "magnetic universe", where the electric compo- 
nent of the electromagnetic field is assumed vanishing. 
A concrete model, previously studied in [ll], [l^] , where 
positive and negative powers of F are included in the 



NLED Lagrangean, is investigated in section IV. In sec- 
tion V we explore the impact of braneworld gravity on the 
dynamics of power-law NLED-magnetic universes. We 
include Randall-Sundrum, as well as Dvali-Gabadadze- 
Porrati brane models [H, Section VI is dedicated 
to the physical discussion of the results obtained, and 
to the statement of brief conclusions. In order to make 
the paper self-contained, an appendix with brief tips on 
how to apply the dynamical systems tools to situations 
of cosmological interest is included. Here we use natural 
units (SttG = Stt/top; = h = c = 1). 



II. NONLINEAR ELECTRODYNAMICS 
COUPLED TO GENERAL RELATIVITY 

The four-dimensional (4D) Einstein-Hilbert action of 
gravity coupled to NLED is given by 



S ^ J d*xy/^[R + L^ + L{F,G)], (1) 

where R is the curvature scalar, L-y - the background 
perfect fluid's Lagrangean, and L(F, G) is the gauge- 
invariant electromagnetic Lagrangean, which is a func- 
tion of the electromagnetic invariants F = F^'^Ff^^i, and 
G = \ea/3fj,uF°'^F^i, [ij]. As usual, the electromagnetic 
tensor is defined as: F^i, — A^^p_ — A^^^u (here the comma 
denotes partial derivative in respect to the spacetime co- 
ordinates, while the semicolon denotes covariant deriva- 
tive instead). Standard (linear) Maxwell electrodynamics 
is given by the Lagrangean L{F) = —F/A. The corre- 
sponding field equations can be derived from the action 
(HI) by performing variations with respect to the space- 
time metric g^^, to obtain: 

Tn-y I rjiEAI 

'^fJ.U — -i/j j/ I J- fill 1 

where 



T'^u = {Pi + Pi) Uf^Ui' - Pi9t^iy, (2) 
T^J"' = 5m- [L{F) - GLg] - AF^aF, " Lf, (3) 

with = P-yit), P-y — P-yit) ~ the energy density and 
barotropic pressure of the background fluid, respectively, 
while Lp = dL/dF, Lpp = d'^L/dF^, etc. Variation 
with respect to the components of the electromagnetic 
potential A^^ yields to the electromagnetic field equations 



F^'' LF + -e''^'"'F^pLG] -0. (4) 
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In this paper we shall consider a homogeneous and 
isotropic Friedmann-Robertson- Walker (FRW) metric of 
the form 



dr^ 



1 — kr'^ 



where a{t) is the cosmological scale factor, and the con- 
stant k = —1,0,-1-1 parametrizes the curvature of the 
spatial sections. In order to meet the requirements of 
homogeneous and isotropic cosmologies (as, in particu- 
lar, the one associated with FRW spacetimes), the en- 
ergy density and the pressure of the NLED field should 
be evaluated by averaging over volume. To do this, we 
define the volumetric spatial average of a quantity X at 
the time t by (for details see [l3] and references therein) : 



III. MAGNETIC UNIVERSE BASED ON 
POWER-LAW NLED 

For the purposes of the present investigation we shall 
consider electromagnetic Lagrangean densities that de- 
pend on the invariant F only. As mentioned in [l^ . a 
particularly interesting case arises when only the aver- 
age of the magnetic part B is different from zero, leading 
to the so called magnetic universe. This case turns out 
to be relevant in cosmology as long as the averaged elec- 
tric field E is screened by the charged primordial plasma, 
while the magnetic field lines are frozen jT5| . 

In this simpler case, the cosmological equations can be 
written in the following form: 



X = 




(5) 



where V = J d^x^—g and Vq is a sufficiently large time- 
dependent three-volume. Following the above averaging 
procedure, for the electromagnetic field to act as a source 
for the FRW model we need to impose that (the Latin 
indexes run over three-space); 



E, = 0, 0, E^Bj = 0, 



EiE-i 



-E gij, BiBj 



(6) 



(7) 



Additionally it has to be assumed that the electric 
and magnetic fields, being random fields, have coher- 
ent lengths that are much shorter than the cosmologi- 
cal horizon scales. Under these assumptions the energy- 
momentum tensor of the electromagnetic (EM) field - 
associated with the Lagrangean density L = L{F, G) - 
can be written in the form of the energy-momentum ten- 
sor for a perfect fluid [12]: 



TnJ^ = {PB +Pb)u^Uu -psg^^u, 



where 



Pb = -L + GLg~4.LfE^, 
Pb = L-GLg-^ (2B2 _ 



E')L, 



(8) 

(9) 
(10) 



E and B being the averaged electric and magnetic fields, 
respectively. In what follows, to simplify the analysis, we 
shall consider a FRW universe with fiat spatial sections 
{k ~ 0), filled with a "magnetic fluid", i. e., the elec- 
tric component E will be assumed vanishing. Even this 
simplified picture can give important physical insights. 



3i/^ = p-y- L, 2H = -jp^ - 
p^ + 3Hjp^ ^0, F + AHF : 



-FLi 



3 
0, 



(11) 



where H = d/a is the Hubble parameter, 7 is the 
barotropic index of the background's perfect fluid {p^ = 

(7 - 1) /O7)' '^liile 
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PB ^ -L , PB ^ L - -^FLf , 

and, as already mentioned, we are considering F ~ 2B^. 
The system of cosmological equations ([Tl]) is a closed 
system: there are 4 unknown field variables {H, B, L, p), 
and 4 differential equations. However, as with any system 
of non-linear (second-order) differential equations, it is 
a very difficult (and perhaps unsuccessful) task to find 
exact solutions. In such cases the dynamical systems 
tools come to our rescue. These very simple tools give 
us the possibility to correlate such important concepts 
like past and future attractors (also saddle equilibrium 
points) in the phase space, with generic solutions to the 
set of equations (|lip . without the need to analytically 
solve them.^ 

In order to be able to apply these tools to the case of 
interest in this paper, we introduce the following phase 
space variables: 



3iJ2 



y = 



FLp 
3IP 



which allow to translate the original problem from the 
space of the field variables {H, B, L, p) into the {x,y)- 
phase plane. After the above choice of variables, the 
system of equations (jlip can be written in the form of 



^ For concise recipes of how to apply these tools to situations of 
cosmological interest we refer the reader to the section Appendix. 
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the following autonomous system of ordinary differential 
equations (ODE): 




-2y 



H 



my 



H 3 , , „ 
— = --7(l-x)+2y 



(12) 



where we have defined the new parameter m — 
LppL/L^, and the prime denotes derivative with respect 
to the new time variable r = In a - properly the number 
of e-foldings. The above autonomous system of ODE is 
not, in general, a closed system of equations, unless the 
parameter m is restricted to be a constant. This will be , 
precisely, the case we will be considering here. After mak- 
ing appropriate physical considerations, the phase space 
corresponding to the above autonomous system of ODE 
can be written as, 

* {(x,2/) \ <x <1, y e R} . 



Since we are considering constant values of the param- 
eter m = LpF ^ 
expression, we obtain: 



L/Lj^, then, after integrating the above 



(13) 



where Lq is an integration constant - which has to be set 
equal to —1/4 in order to be able to recover the standard 
Lagrangean for Maxwell ED when m — -, and 



L = Lq exp (XF), 



1, 



where Lq and A are integration constants. 

For the cases of interest m — const the phase space 
variables x, y are not independent variables since their 
ratio is a constant: 



X 

y 



FLf 



1 . 



Hence the phase space is just a one-dimensional space, 
and the autonomous system of ODE reduces to just 
one independent equation (m ^ 1): 



3(171- 1)7-1-4 

TO — 1 



x{x — 1). 



(14) 



The relevant magnitudes of observational interest can 
be given in terms of the variable x and of the parameter 
to: 



flB = X, fl-y — 1 ~ X, ojb — 



3m -I- 1 
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3(1 - to) 
3(m- 1)7 + 4" 
2(to- 1) 



Note that the equation of state (EOS) parameter of the 
NLED component a;^, is the same value no matter what 
the critical point is. It depends only on the value of the 
overall constant to. For to = the NLED fluid behaves 
just like standard electromagnetic radiation with lob = 
1/3 (in fact linear electrodynamics), while as to 00, 
ujb —1: the NLED fluid behaves as a cosmological 
constant (quantum vacuum). 

The equilibrium points corresponding to this case are 
obtained by requiring that in Eq. (1141) a;' = 0. Two 
equilibrium points are found: 

• The matter (background fluid)-dominated solution 
{x = ilB = 0): 

37 

n^ = i, g = -i + _i. 

• The NLED-dominated solution (x — ^Ib = 1): 

il-y — U, q — —- 



TO - 



1 



Note that, for the case to = 1 (exponential Lagrangean 
above), since x/y — always, only the matter-dominated 
equilibrium point is found. In this case the EOS param- 
eter of the magnetic fluid is undefined. This case will not 
be considered anymore in this paper. 

To explore the linear stability of these critical points 
X = Xc, we expand the equation (|14l) in their neighbor- 
hood and keep terms linear in the small perturbation e. 
For the first equilibrium point x ~ Xc — 0, we obtain the 
following equation for the perturbation: 



e' = fee, k — 



3(to- 1)7-^4 



1 



while, for the second critical point (x — x^ — 1): 

e' ^-fc(e + l) . 

We see that, for positive fc > the matter-dominated 
solution is unstable, while the NLED-dominated phase is 
stable. For fc < these equilibrium points interchange 
their stability properties: the NLED-dominated phase is 
unstable while the matter-dominated solution is stable. 
Consider, for instance, background dust with 7 = 1. In 
this case 

i + l 



to — 1 
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so that, for TO > 1 the matter-dominated solution is 
the past attractor, while the NLED inflationary solu- 
tion {ujb < 0, g < 0) is the late-time attractor. For 
TO < 1, on the contrary, the NLED non-inflationary so- 
lution {ujb > 0, q > 0) is the past attractor while the 
matter-dominated - also non-inflationary {q — 1/2) - so- 
lution is the future attractor. 

One then expects that a combination of NLED terms 
proportional to positive and negative powers of F, would 
eventually allow for a cosmological scenario where the 
universe might evolve from a past NLED-driven (non- 
inflationary) state into a future (late-time), also NLED- 
driven, inflationary stage, transiting through a matter- 
dominated solution. The study of such a potentially in- 
teresting scenario will be the subject of the next section. 



IV. UNIFIED DESCRIPTION OF BOUNCING 
COSMOLOGY FOLLOWED BY LATE-TIME 
ACCELERATED EXPANSION 

In this section we focus in the study of the asymptotic 
properties of a cosmological model with interesting fea- 
tures, namely a phase of current cosmic acceleration and 
the absence of an initial singularity, which was proposed 
in [ll] (see also [ll|) and is based upon the following 



Lagrangean density: 



L = — F 
4 



aF^ + PF- 



(15) 



where a and /3 are arbitrary (constant) parameters. As 
seen this Lagranagian contains both positive and nega- 
tive powers of F. The second (quadratic) term dominates 
during very early epochs of the cosmic dynamics, while 
the Maxwell term (first term above) dominates in the ra- 
diation era. The last term in (jTS]) is responsible for the 
accelerated phase of the cosmic evolution fT3| . The above 
Lagrangean density yields a unified scenario to describe 
both the acceleration of the universe (for weak fields) and 
the avoidance of the initial singularity, as a consequence 
of its properties in the strong-field regime. 

Recalling that we are considering magnetic universes, 
i.e., F = 2B^, where is an averaged value of the mag- 
netic field, ^ the stress-energy tensor associated with ([15]) 
can be written in the form of an equivalent perfect fluid 
stress-energy tensor with energy density and parametric 
pressure (see section II): 



PB = -(l-8aB^-A 



(16) 



respectively. Notice that, for large values of the NLED 
field, positivity of energy requires that B < l/-\/8a, 
while, for small enough values of i? <C 1, if one consid- 
ers positive /3 > 0, then positivity of energy implies that 
B > {7py/^. In the latter case, in the unified model 
given by the Lagrangean density ([T5]) [T2j . there arise 
both higher and lower bounds on the values the NLED 
field B can take. The existence of the lower bound, at 
first sight might appear problematic, however, given that 
the observational data constraints the parameter /3 to be 
« 4 X 10~^* g cm~^ [Tl|, then the lower bound 
on B can be admitted without going into conflicts with 
observations. 

In the rest of this section we assume the background 
fluid to be dust cold dark matter (CDM), so that 7 = 1. 
Our goal in this section will be to put the corresponding 
cosmological equations dTTI) : 



= Pcdm + Pb , 
= -pcdm - {pB +Pb), 
Pcdm + '^Hpcdrn ~ 0, 

PB + 3H{pB +Pb) = => B = -2HB, (17) 



with Pb and pb given by (|16p . in the form of an au- 
tonomous system of ODE. In the present case, since there 
exists a combination of powers of F in the Lagrangean, 
we can not use the variables of the former section. We 
choose the following phase space variables instead: 



PB WaB" 
X=-rT^, y = , z 



4/3 



In terms of these variables one has: 

, y z pB+PB . 

X , = 4:X — V + Z , 

4 8' '^^ ' 



(18) 



and, also 



2^ = -3-x + y-z, 

meanwhile, the dimensionless energy density of the CDM 
can be expressed as function of the x-variable alone 
i^B^x); 



^ For details of the averaging procedure consult [12 
section II of the present paper. 



See also 



j-j Pcdm -| 

^^cdm — ^j^2 



(19) 
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Other magnitudes of observational relevance are the 
equation of state (EOS) parameter of the magnetic field, 
and the deceleration parameter: 

X — y + z 1 + X — y + z 



respectively. 

The following autonomous system of ODE is obtained 
out of ([T71): 

x' — {x — l){x — y + z), 
y' = -y(5 ~x + y- z), 
z' = z{7 + x-y + z). (20) 

Here, for generality of the analysis we shall consider 
arbitrary a G and /? € SR. The phase space relevant to 
the present study is then given by the following region in 
{x,y,z): 

'ifu = {{x,y,z)\Q<x<l, (y,z)e5R2, 

8x + 2y + z>0}, (21) 

where we have taken into consideration that 
B'^ _ y z 

The equilibrium points of pO)) in ^c/, their existence 
and stability properties, are discussed below: 

1. Radiation-dominated phase: 

Prad={x,y,z) = (1,0,0), Hb^I, n,dm = 0. 

This is a decelerating expansion solution {q — 1), 
that is fueled by standard radiation with lob — 1/3- 
The eigenvalues of the linearization matrix corre- 
sponding to this equilibrium point are: 

Ai = 1, A2 = -4, A3 = 8 , 

so that it is a saddle in '^jj. 

2. Infra-red NLED-dominated solution: 

Pi^^rf- (1,0,-8), ne^l, f^c*n-0. 

This is a late-time, super-inflationary solution (q = 
—3), where the NLED fluid mimics phantom behav- 
ior (lob — —7/3). This solution exist only for z < 
(negative /3 < 0). In this case there is no lower 
bound on the magnitude of the magnetic field. The 
eigenvalues of the linearization matrix for P^^^ are 

Ai = -8, A2 = -7, A3 - -12 , 



so that this solution corresponds to a late-time (fu- 
ture) attractor. Since, this point is dominated by 
the NLED magnetic field mimicking phantom field, 
then, the late-time attractor could be associated 
with a cosmological singularity (most probably a 
big-rip type of singularity). 

3. Ultra-violet NLED-dominated phase: 

This solution corresponds to an early-time, super- 
stiff- fluid solution {ujb — 5/3), which is associated 
with super-decelerating expansion (q — 3). This 
solution exist only for y < 0, i. e., if the con- 
stant a is a negative quantity (a < 0). Curi- 
ously, this case does not meet the conditions for 
a bounce (there is no upper bound on the magni- 
tude of the magnetic field), and the corresponding 
cosmology starts with a big-bang singularity, since 
B^/H^ ^ 0, B ^ 0, ^ H ^ 00. This solution is 
a past attractor in ^1^^;, since the eigenvalues of the 
corresponding linearization matrix: 

Ai = 5, A2 - 12, A3 = 4, 

are all positive quantities. This means that this 
point is the starting point of every probe path in 
the phase space ^[/. 

4. CDM-dominated solution: 

Pcd,n = (0,0,0), CtB^O, rtcdrn^l. 

This phase of the cosmic evolution is characterized 
by decelerated expansion (q — 1/2). The EOS pa- 
rameter for the magnetic field is undefined in this 
case. The existence of this solution is necessary for 
the formation of the observed amount of structure. 
This is also a saddle equilibrium point in since: 

Ai = -5, A2 - 7, A3 - -1 . 

For positive definite a and (3, only the radiation- 
dominated equilibrium point Pnied, and the CDM- 
dominated solution Pcrfm, are found in ^£/. 

The above results confirm our expectation that a com- 
bination of positive and negative powers of the electro- 
magnetic invariant F in the NLED Lagrangean, can drive 
a very interesting cosmological scenario leading to accel- 
erated expansion at late times. Unfortunately, early-time 
(primordial) inflation can not be obtained in the present 
model. It is remarkable also that the bounce is not a 
generic solution in such models. 

In the next section we shall explore the influence of 
the brane effects on the dynamics of the NLED-magnetic 
universe. 
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V. NLED ON THE BRANE 

In this section we explore the possible effect of 
braneworld gravity on the above picture. A Randall- 
Sundrum braneworld model Ts*] has proved satisfac- 
tory to modify the primordial inflation scenario, since 
the extra-dimensional brane effects are appreciable at 
high-energies (large brane tension) . A Dvali-Gabadadze- 
Porrati braneworld model, on the contrary, modifies 
gravity at large (cosmological) scales [l^l- Hence, 
the study of the cosmological dynamics driven by 
NLED-magnetic universes trapped in either RS or DGP 
braneworlds, is an interesting possibility to look for al- 
ternative, viable, cosmological scenarios. 



A. RS2 brane 

One of the most appealing braneworld scenarios is the 
Randall-Sundrum brane model of type 2 (RS2) (12]. In 
this model a single co-dimension 1 brane with positive 
tension is embedded in a five-dimensional anti-de Sitter 
(AdS) spacetime, which is infinite in the direction per- 
pendicular to the brane. In general, the standard model 
particles are confined to the brane, meanwhile gravita- 
tion can propagate in the bulk. In the low-energy limit, 
due to the curvature of the bulk, the graviton is confined 
to the brane, and standard (four-dimensional) general 
relativity laws are recovered. 

The FRW cosmological equations for a RS2 brane with 
CDM and a perfect fluid of NLED trapped on it, can be 
written as; 



2H 



Pcdrn 



- Pcdm + PB, 



(22) 



where A is the brane tension. 

In order to write the above system of equations into the 
form of an autonomous system of ODE, we introduce the 
following phase space variables (recall that ps = —L); 



-L 



r, y 



After this we can write, 

PT 2(1 -t-) 



Pt 



(23) 



A 



Notice that, in the above phase space variables, the high- 
energy regime of the brane is associated with the limit 



f — > 0, while the low-energy, general relativity limit cor- 
responds to t; = 1. 

The following autonomous system of ODE is found: 



x' = 2\ 2y~x 



H 



v' ^ -3{v~x) + ^y~2v^, (24) 



where 



H 



3(2- w) 



4y] , 



i/2 2v ' 

and, as before, for m ^ 1 (see section III), 

m = - const. ^ L = LoF^d-™) . 

Besides, since y = x/{m—l), only one of these variables is 
independent, say, x. This means that, in fact, the system 
of ODE ((M)) is a two-dimensional system of ODE: 




3(to- 1) 



2v 



H 



(25) 



with 



^ H 



3(2-1;) 



m — 1 



The phase space where to look for the equilibrium points 
of the above autonomous system of ODE is given by the 
following compact {x, t;)-region: 

*6rane = {{x,v)\0 <X<1, Q<V<1). 



The equilibrium points of (|25|) . their existence and sta- 
bility properties are discussed below. Only two equilib- 
rium points can be found: 

• CDM-dominatcd (GR) critical point: 

P^d^ = {x,v)^{{),l), r!eA„ = l, f^s=0. 

This point is associated with decelerated expansion 
q = 1/2. The eigenvalues of the corresponding lin- 
earization matrix are: 



Ai 



-3, A2 



3rn -|- 1 
m — 1 
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For 



1. Combination of positive and negative powers of F 



-- < TO < 1 , 

the point P^^^ G "^brane is the late-time attractor, 
while for 

TO < — - , or, TO > 1 , 



• Scaling (RS) critical point: 

P^c^ = (3fcuo, Wo), ^^cdm = (1 - 3fc)uo, = 3fcuo 

where 

'3to + 1\ 57to+11 



k = 



9m — 5 



3(llm+ 1) ■ 

The scaling character is due to the existence of non- 
vanishing ratio 

^cdm _ I -3k _ 8 

~ 3fc ^ 9TO-I-3 ■ 

This solution exists only for m < —1/3. The decel- 
eration parameter is given by: 

TO + 1 

q = 7 • 

TO — 1 

The eigenvalues of the Jacobian matrix correspond- 
ing to this critical point are: 

3to+1 4(57TO-hll) 

Al — A2 — 



TO — 1 ' ^ (9to — 5)(to — 1) 

It can be verified that, for the TO-interval where the 
point exists, both eigenvalues are negative quanti- 
ties, so that the scaling solution is always the late- 
time attractor. 

Note that for to > —1/3, only one equilibrium point 
(the CDM-dominated solution) exists, while for to < 
— 1/3 both critical points exist. In the latter case, the 
CDM-dominated solution is a saddle critical point in 
"^brane, wliilc the Scaling solution is the late-time attrac- 
tor. The latter attractor is associated with accelerated 
expansion whenever to < — 1, while for — 1 < to < —1/3, 
the expansion is decelerated instead. Notice, also, that 
the NLED field EOS parameter ujb ■= -1 - 4/3(to - 1), 
does not depend on x, z, so that is the same value for 
both the CDM-dominated solution and for the scaling 
one. 

If we compare the above results with the ones of sec- 
tion III, we see that the RS2 brane does really modify 
the dynamics of the NLED-magnetic universe. Actually, 
the NLED-dominated solution reported in section III has 
been replaced by the scaling solution above. Hence, it 
makes sense to inspect the model of section IV when the 
brane effects are included. 



Here we shall consider the Lagrangean ([T5)) which con- 
sists of a combination of positive and negative powers of 
the electromagnetic invariant F. The cosmological equa- 
tions are the equations (j22p , which can be rewritten in a 
less formal way: 



Pcdm + PB, 



2H 



-{Pcdm + PB +Pb) (1 + 



PT 
A 



Pcdm — —iHpcdm, B — -~2HB, 



(26) 



where pB and pB are defined in equations (jl6p . It is con- 
venient to introduce the following phase space variables: 

_ PB _ l6aB^ 4/3 _ PT 

^-3i/2' y- 3i/2 ' ^-3^/252' ^-3i/2- (^') 

In terms of these variables one has, for instance: 

'1 - v' 



PT _^ 

A " 



and, also, 



H' , , f2-v 

2— = -{3v + x-y + z) 

H \ V 



The variable v controls the brane regime, so that, for 
instance, the GR-limit (formal limit A — >■ 00) corresponds 
to w = 1. The following autonomous system of ODE can 
be derived out of (l26l): 



(x-y + z) 
{x~y + z) 
{x~y + z) 



X = 3(1 — v)x + 

-5 + (l-t; 



7+il-v) + 



w' = (1 — f)(3w + X — y + z). 



(28) 



At f = 1 the first three equations above coincide with 
the equations (|20|) of section IV, which hold for general 
relativity with a NLED-magnetic field. The phase space 
of the model can be defined as follows: 



^U-^- = {{x,y,z,v)\0<x<l, (2/,z)e5R2, 

8a; + 2y + z>0, < w < 1}. (29) 
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For the following discussion it will useful to write sev- 
eral magnitudes of observational interest in terms of the 
above phase space variables = x): 



LOB 



X — y + z 
3x 



4- 3w 



2v 



(x-y + z) 



The critical points of the autonomous system of ODE 
in the phase space vj/^""^^ their physical proper- 
ties and stability are discussed below. As in the general 
relativity case (section IV), four equilibrium points are 
found: 



1. CDM-dominated solution: 



Pr, 



(0,0,0,1), r!s = o, n^dm^i 



Since q = 1/2, this solution is associated with decel- 
erated expansion. The NLED-magnetic field EOS 
parameter wb is undefined. This is a saddle equilib- 
rium point in x]/^''"'=. Actually, the eigenvalues of 
the linearization matrix corresponding to this point 
are: 

Ai = —5, A2 = 7, A3 = — 3, A4 = —1 . 



2. Radiation-dominated solution: 

^rad = (1, 0, 0, 1), fis = 1, r^cdm = . 

It is also a decelerated-expansion solution {q = 1) 
driven by standard (Maxwell) radiation {ub = 
1/3). As the CDM-dominated phase, this solution 
also represents a saddle equilibrium point in \I'^°"'^ , 
since the eigenvalues of the linearization matrix are 
of opposite signs: 



Ai — 8, A2 



1, A 



3,4 



3. UV NLED-dominated solution: 



Cied = (l>-4,0,1), Ob = 1, n,am = 0. 

This solution shares many properties with its sim- 
ilar in section IV: it is a super-decelerated (g = 3), 
super-stiff state (ujb = 5/3), associated with a big- 
bang- type singularity {H -^00). The new feature 
here is that the stability properties have been mod- 
ified by the brane effects. Actually, the eigenvalues 
of the linearization matrix in the present case are: 

Ai ^ 12, A2 =4, A3 = 5, A4 = -8 , 

so that this is a saddle equilibrium point in vl/^«"<= 
(its similar in section IV is a past attractor). 



4. IR NLED-dominated solution: 



-Pnied^ (1,0,-8,1), nB = l, n,dm=0. 

As in the former case, this solution shares many 
properties with its similar in section IV: it is a 
super-accelerated (g = —3), phantom-like solution 
{ujb = —7/3), possibly associated with a big- rip- 
type singularity. Other properties, the stability in 
particular, have been modified by the brane effects. 
Actually, the eigenvalues of the Jacobian matrix are 

Ai = -12, A2 = -8, A3 = -7, A4 = 4 , 

so that it is also a saddle equilibrium point in the 
phase space (PU)) . 

The first thing that is worthy of mention, is the fact 
that all of the above critical points represent saddle equi- 
librium points in the phase space v]>^<^"'= (j29p . Note 
that, since v = 1, the four equilibrium points are as- 
sociated with GR. In fact these coincide in almost all 
aspects with the ones found in section IV. The first two 
points Pcdm, and Prad, show no fundamental differences 
with their similar in that section. However, the NLED- 
dominated solutions P^^d ^^"^ ^nied^ have different sta- 
bility properties than their corresponding solutions for 
the GR case: while in the latter case P^ied ^as the 
past attractor and P^i^d future attractor, in the 

present case both are saddle equilibrium points as already 
said. This means, in turn, that the space-time singular- 
ities associated with these critical points (big-bang and 
big-rip singularities respectively), might be evaded in the 
present case. Modification of the stability properties of 
the UV solution was expected since RS brane effects are 
appreciable at early times (high energies). The surprise 
was the IR solution since, it is expected that at low ener- 
gies/large cosmological scales, the RS2 brane effects can 
be safely ignored. We postpone the discussion of the 
possible origin of this effect to section VI. 



B. DGP brane 

Another braneworld scenario that has received much 
attention in the last years is the Dvali-Gabadadze-Porrati 
(DGP) model [l3l- It describes a brane with four- 
dimensional (4D) world-volume, that is embedded into 
a flat 5D bulk, and allows for infrared (IR)/large-scale 
modifications of gravitational laws. A distinctive ingredi- 
ent of the model is the induced Einstein-Hilbert action on 
the brane, that is responsible for the recovery of 4D Ein- 
stein gravity at moderate scales, even if the mechanism 
of this recovery is rather nontrivial [T^ . The acceleration 
of the expansion at late times is explained here as a con- 
sequence of the leakage of gravity into the bulk at large 
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(cosmological) scales, so it is just a 5D geometrical ef- 
fect. As with many IR modifications of gravity there are 
ghosts modes in the spectrmxi of the theory [l3l-'^ Never- 
theless, studying the dynamics of DGP models continues 
to be a very attractive subject of research. 

The cosmological FRW equations for a DGP brane 
with CDM and NLED-magnetic fluid trapped on it are 
the following: 



value V — 1. This subspace corresponds to general rela- 
tivity. For the self-accelerating branch (minus sign in the 
definition of Q) the phase space is given by 

*dsp = {(^-")|0<^<l, 0<t;<l}, (33) 

while, for the Minkowski cosmological phase (plus sign in 
the definition of Q) it is given by 



iQl = {Pcdrn + Pb) , Q'i = H^± 

Pm^-3Hp,d^, F + 4HF = 0. 



H 



(30) 



In order to write ([50)) as an autonomous system of ODE 
we introduce the following phase space variables: 



L 



Q 
H' 



(31) 



The following system of ODE is found: 



/3to+ 1 
V 1 — m 



{x — l)x, 



V = 



1 



o /3m -hi 

3-h I X 

\ 1 — n 



(32) 



In terms of x, v, the relevant parameters can be written 
as: 



cdm 



v'^{l — x), V,B—V^X, UJb = 



and the deceleration parameter 



2 

1 - - 
3 



{3m + l)x-3{m- 1) 



m — 1 



m + 1 

1—771 ' 



1 + 



Note that the EOS parameter of the NLED-magnetic 
fluid Ldg depends only on m and not on the concrete 
critical point being considered. For 777 > 1 the NLED 
component mimics phantom fluid {ljb < ^ 1)- 

According to which branch of the DGP is being con- 
sidered, the phase space for this case is split into two 
regions that share a common subspace deflned by the 



^ In fact there are ghosts only in one of the branches of the DGP 
model; the so-called self-accelerating branch, or self-accelerating 
cosmological phase [T3l . The Minkowski cosmological phase is 
free of ghosts. 



dgp 



{{x,v)\Q <X<1, V>1}. 



(34) 



Points with negative w-s (u < 0), are not considered due 
to the symmetry of the field equations (|30l) under the 
change Q — >■ — Q. 

There are found two equilibrium points belonging in 
the common GR-boundary ^'^gp ^1 ^dgp- 

• CDM-dominatcd solution: 



jGR 



(0,1), r!cdm = l, Ob-0, q- 



1 



Since the eigenvalues of the Jacobian matrix are 

3 3777-1-1 

Ai — -, A2 — — , 

2 777 —1 

then, this point is a saddle (A2 < 0) for 777 < —1/3 
or 777 > 1. For — 1/3< 777<1, itisa source critical 
point/past attractor in ^'^gp ^ ^dgp instead. 

• NLED-dominated solution: 



P„^f, = (l,l), l^cdrn^O, 17b = 1, = - 



777 + 1 
777 —1 



This solution is inflationary {q < 0) whenever |7r7| > 
1. It is the late-time attractor if 777. > 1, since 
the eigenvalues of the corresponding linearization 
matrix are: 

3777+1 2 

Ai — — , A2 — — - 



777 — 1 ' " 777 —1 

The NLED-dominated solution can be also a past 
attractor if —1/3 < 777 < 1. Otherwise it is a saddle 
point in VI/- ^ n VI/+ ^. 

while a similar couple of equilibrium points are found 
in ^'^gp - the phase space corresponding to the self- 
accelerating DGP brane: 

• Self-accelerating solution: 



(0,0), ricdm-0, fis = o, q = -i. 



self — a 

which linearization matrix has the following eigen- 
values: 

3 3777 +1 

Ai — -7:, A2 — — . 

2 777 —1 
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• Alternative self-accelerating solution: 

^'ei/-a = (1,0), f^crfm^O, l^S = 0, ^ = -1, 

whose Jacobian matrix has the eigenvalues: 

3m + 1 2 
Ai — , A2 — . 

m — 1 m — 1 

Both self-accelerating solutions are inflationary and, 
since v — 0, they correspond to the late-time de Sit- 
ter phase, when Q = H = l/vc- The first self- 
accelerating equilibrium point P^^ij:_^ is the late-time 
attractor if— 1/3 < m < 1. Otherwise (m < —1/3, 
m > 1) it is a saddle equilibrium point in "^dgp- Alterna- 
tively, the second self-accelerating critical point Pg^if^^i 
can be either, i) a saddle equilibrium point if m > —1/3, 
or ii) the future attractor if m < —1/3. This means 
that, for m < 1, there always coexist two de Sitter (self- 
accelerating) solutions, one being a saddle equilibrium 
point, while the other one is the late-time (future) at- 
tractor. For TO > 1 the self-accelerating solution is a 
saddle critical point. 



As before 

and, additionally: 

'?=-1+(y^) (3 + x-y + z), 

and now the NLED-field EOS parameter does indeed de- 
pend on the equilibrium point under consideration: 

X — y + z 

According to whether branch of the DGP model one is 
being considering, the phase space can be defined as 

8.x + 2y + z > 0, < t; < 1}, (38) 
for the self-accelerating branch, or 



1. Combination of positive and negative powers of F 

The next step is explore the effects of the DGP brane 
by considering the concrete Lagrangean (|15p . which con- 
sists of positive and negative powers of the electromag- 
netic invariant F. The FRW field equations are: 

3Q^ - {Pcdm + Pb), Q\ = H^±- 

^ = -^-ff {Pcdm + PB+ Pb) , 

p^ = -iHp,d,n. B + 2HB = 0, (35) 

where, as before, we are considering CDM and NLED- 
magnetic fluids on the brane, and pb, Pb are given by 
equations ([TBI) . If introduce the phase space variables: 

_ PB _ l6aB* 4/3 _Q 

''^3Q2' y- 3Q2 ' ^-3g252' ^''^> 

then, the field equations (j35p can be written in the form 
of the following system of ODE: 

x' — {x — l){x — y + z), 
y' = 2/(-5 + x-y + z), 
z' = z{7 + X — y + z), 

.' = -|(i^)(3 + x-y + z). (37) 



8x + 2y + z > 0, w > 1}, (39) 

for the Minkowski cosmological phase instead. As before, 
due to the symmetry of the field equations (I35p under 
the replacement Q — > — Q, critical points with negative 
values of the v-s {v <0), are not being considered. 

The equilibrium points of (j37p in the interception sub- 
set ^ *<lgp (GR-boundary), are: 

1. CDM-dominated solution: 

P,<^^ = (0,0,0,1), n^dm^l, f^B = 0, g=^, 

with the following eigenvalues of the Jacobian ma- 
trix Ai = 3/2, A2 — 7, A3 = —1, and A4 — —5. 
Hence, as before, this point is a saddle in the GR- 
boundary. 

2. Radiation-dominated phase: 

p,^f = (1,0,0,1), r!c*n = o, ris^i, q = i. 

The radiation-dominated solution {ujb ~ 1/3), is 
also a saddle in the phase space, since the eigen- 
values of the linearization matrix are of opposite 
signature: Ai = 1, A2 = 2, A3 = 8, A4 = —4. 

3. UV NLED-dominated solution: 

^'n'^e^'d = (1,-4,0,1), f^cd™ = 0, 17b = 1, g-3. 
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Here, as in the GR case, the super-stiff {lob — 5/3), 
decelerating solution, is the past attractor, since 
Ai = A2 = 4, A3 = 5, A4 = 12. It is associated 
with a big-bang cosmological singularity {B^/H^ = 
^ H ^ 00). 

4. IR NLED-dominated phase: 

P,(,^rf = (1,0,-8,1), n,dr,^^0, ne^l, 9 --3. 

This is the super- accelerating (phantom-like) solu- 
tion, possibly associated with a big-rip-type singu- 
larity, as in the standard GR case. As in the pure 
GR situation, this equilibrium point is the late-time 
attractor since: Ai = —2, A2 = —7, A3 = —8, 
A4 = -12. 

Meanwhile, in the phase space ^^g„, corresponding to 
the self-accelerating branch of the DGP, one founds the 
following critical points: 

1. Self- accelerating phase 1: 

nW-a = (0,0,0,0), n,dm = nB = o, q = -i. 

The eigenvalues are: 
Ai = 7, A2 = -3/2, A3 - -1, A4 - -5 . 

2. Self- accelerating phase 2: 

^.Va = (I'O'O'O)' f^cdm-f^s=0, g = -l, 
with the following eigenvalues: 

Ai = 1, A2 =8, A3 = -2, A4 = -4. 

3. Self- accelerating phase 3: 

^iLi/-a = (1,0,-8,0), r!e*„ = r!B = o, q = -i. 

Here, 

Ai =2, A2 = -7, A3 = -8, A4 - -12 . 

4. Self-accelerating phase 4: 

= (l'-4,0,0), f7eA„ = r!B = 0, g = -l. 
The eigenvalues of the linearization matrix are: 
Ai=4, A2 = 5, A3 = 12, A4 = -4. 

These four critical points correspond to the inflationary, 
late-time self- accelerating solution H = 1/rc (y = Q 
Q = 0), and are saddle equilibrium points in ^^np- 

As seen, the Minkowski cosmological phase of tne DGP 
does not modify the stability properties of the GR equi- 
librium points found in section IV. The self-accelerating 
branch of the DGP differs from results in that section in 
the addition of four critical points, associated with the 
self-accelerating de Sitter solution H = 1/rc- 



VI. DISCUSSION AND CONCLUSIONS 

Non-linear electrodynamics can supply a useful sce- 
nario where to discuss such relevant problems of the stan- 
dard cosmological model as the initial singularity related 
with the big-bang, and the early-times and present stages 
of accelerated expansion of the universe. A Lagrangean 
density of the form given in (jlSp . for instance, can ac- 
count for the present accelerated stage of the cosmic evo- 
lution, as well as for avoidance of the initial singularity, 
in a unified picture p^ . In a general setting, where arbi- 
trary powers of F are invoked, this behavior can be ex- 
plained as due to the asymptotic properties inherent in 
these kinds of Lagrangean. Actually, given a Lagrangean 
of the form 

L = LoF^/^^~"''> , 

it has been shown in section HI, that, for m > 1, the 
matter-dominated solution is the past attractor, while 
the NLED inflationary solution {ujb < 0, g < 0) is 
the future/late-time attractor. For m < 1, these so- 
lutions exchange their stability properties: the NLED 
non-inflationary solution (ws > 0, g > 0) is the past at- 
tractor, while the matter-dominated solution (also non- 
inflationary since q — 1/2), is the late-time attractor. 
While none of these models alone seems like a plausible 
cosmological model, combinations of negative and posi- 
tive powers of the electromagnetic invariant F, provide 
an interesting cosmological scenario instead. 

To make the analysis simpler, we have found instruc- 
tive to consider the so called "magnetic universe" (no 
electric component). For deflniteness we have explored 
a FRW universe whose dynamics is based on the La- 
grangean 

Ltot = Lg + Lcdrn + L{F) , 

where Lg is the standard Einstein-Hilbert Lagrangean, 
Lcdrn corresponds to a pressureless perfect fluid (CDM), 
and L{F) is the NLED Lagrangean displayed in Eq. (fTS)) . 
By means of the dynamical systems tools we are able to 
find the relevant asymptotic dynamics of the model. The 
goal is to write the original system of cosmological equa- 
tions in the form of an autonomous system of ordinary 
differential equations. Then, the information on the dy- 
namics is encoded in the structure of the phase space, 
in particular in the knowledge of the equilibrium/critical 
points and their stability properties. The knowledge of 
the equilibrium points in the phase space corresponding 
to a given cosmological model is a very important infor- 
mation since, independent on the initial conditions cho- 
sen, the orbits of the corresponding autonomous system 
of ODE will always evolve for some time in the neigh- 
borhood of these points. Besides, if the point were a 
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stable attractor, independent on the initial conditions, 
the orbits will always be attracted towards this point (ei- 
ther into the past or into the future). Going back to the 
original cosmological model, the existence of the equilib- 
rium points can be correlated with generic cosmological 
solutions that might really deside the fate and/or the 
origin of the cosmic evolution. In a sense the knowl- 
edge of the asymptotic properties of a given cosmological 
model is more relevant than the knowledge of a partic- 
ular analytic solution of the corresponding cosmological 
equations. While in the later case one might evolve the 
model from given initial data giving a concrete picture 
that can be tested against existing observational data, 
the knowledge of the asymptotic properties of the model 
gives the power to realize which will be the generic be- 
haviour of the model without solving the Einstein's field 
equations. In the dynamical systems language, for in- 
stance, a given particular solution of the Einstein's equa- 
tions is just a single point in the phase space. Hence, 
phase space orbits show the way the model drives the 
cosmological evolution from one particular solution into 
another one. Equilibrium points in the phase space will 
correspond to solutions of the cosmological (Einstein's) 
equations that, in a sense, are preferred by the model, i. 
e., are generic. The lack of equilibrium points that could 
be correlated with a given analytic solution of the model, 
amounts to say that this solution is not quite generic, oth- 
erwise unstable in terms of phase space variables, and can 
not be taken too seriously. 

The asymptotic structure of the FRW cosmological 
model based on the above Lagrangean Ltot shows that 
the UV NLED-dominated, super-stiff-fluid (decelerating 
expansion) solution, is the starting point of the cosmic 
evolution (all orbits in the phase space emerge from this 
solution). The radiation and CDM-dominated solutions 
are also critical points of the corresponding system of 
ODE, so that there is room in this model to explain 
the observed amount of structure we see. Besides, the 
IR NLED-dominated, super-inflationary solution is the 
future/late-time attractor, meaning that it is the end- 
point of the cosmic evolution. Due to the term cx I3B~^ in 
the NEED Lagrangean the NEED- fluid at late times 
behaves like phantom matter: the IR NLED-dominated 
future attractor could be associated with a big-rip type 
of cosmological singularity. Hence, in the above model 
both, the starting singular point of the cosmic evolution - 
the bib bang singularity -, as well as the catastrophic fate 
of the cosmic evolution, are unavoidable. The fact that 
the past attractor is associated with the initial cosmo- 
logical singularity means that, contrary to what former 
authors found [1, [l3| , avoidance of the big bang singular- 
ity is not a generic property of models of the kind given 
by the Lagrangean (|15p . Only under very specific initial 
conditions, or under very special considerations, solutions 



were the singular origin of the evolution is evaded, may 
be found. 

When the brane effects are taken into consideration, 
one finds that while the DGP brane does not apprecia- 
bly affects the cosmic dynamics driven by the NEED 
Lagrangean (|T5t - but for the occurrence of critical 
points associated with the de Sitter solution H = l/rj, 
in the self- accelerating branch of the DGP model -, the 
Randall-Sundrum brane effects do modify the stability 
properties of the GR equilibrium points. The interest- 
ing finding is that the extra-dimensional effects modify 
not only the UV-asymptotics, but also the IR/late-times 
dynamics. Actually, once the RS2 brane effects are con- 
sidered, the attractor structure of the UV (past attrac- 
tor) and IR (future attractor) NLED-dominated solu- 
tions, within GR, is demoted to just saddle property. In 
other words, the consequence of the extra-dimensional 
(RS2) brane effects is that, both NLED-dominated criti- 
cal points above, represent just saddle equilibrium points 
in the phase space. This would mean, in turn, that the 
big-bang (initial) and the big-rip (final) cosmological sin- 
gularities, can be evaded just by a proper choice of the 
initial conditions: there are orbits of the autonomous sys- 
tem of ODE that do not pass through the neighborhood 
of the corresponding critical points. 

The finding that the stability properties of the infrared 
(IR) NLED-dominated solution are modified by the RS2 
brane effects, came as a surprise. Actually, only at very 
high energies can the graviton acquire large momenta 
along the extra dimension and may escape into the bulk 
(5D) spacetime. On the contrary, at low energies (large 
cosmological scales), it is expected that the RS2 brane 
effects can be safely ignored. However, while making such 
statements one has to be careful. In the cosmological 
context, the most evident RS2 brane effect is to modify 
the Friedmann equation: 

3H».,.(l + g). 

Hence, at very high energy density 3> A (much big- 
ger than the brane tension), the Friedmann equation is 
fundamentally modified 3iJ^ oc p^. If in the course of 
the cosmic expansion the total energy content of the 
universe dilutes, then as long as pT ^ A, one recov- 
ers standard GR-Friedmann behavior. Now look at the 
Lagrangean density for the NLED-magnetic field (fT5|) 
considered here. Note that as the expansion proceeds 
the magnetic field F (x dilutes, and the component 
cx f3B~^ in (IT5|) grows without limit. This means that 
the total energy content of the universe starts growing 
at the expense of the NEED component so that, at late 
times, eventually, pT might become much larger than 
the brane tension once again, rendering the brane ef- 
fects important at late times also. This effect has been 
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demonstrated to be generic of phantom fields of any na- 
ture which are trapped in a RS2 braneworld |l^. Here 
we have demonstrated that the above mixing of UV/IR 
scales, due to extra-dimensional effects on the dynamics 
of phantom matter, is distinctive only of models where 
the brane effects modify the right-hand-side (energy den- 
sity) of the Friedmann equation, so that, for instance, the 
DGP brane model does not show this mixing. 

Summarizing: we have investigated the asymptotic 
structure of magnetic-NLED models where combinations 
of positive and negative powers of the electromagnetic 
invariant F are considered. It has been demonstrated 
that such models can supply an interesting cosmological 
scenario where the end-point of the cosmic evolution is 
a super-accelerated, singular state indistinguishable from 
the big-rip inherent in (scalar) phantom driven cosmolo- 
gies. RS brane effects may drastically change the nature 
of the starting point, as well as the fate of the cosmic 
evolution: both, the big bang singularity, and the big 
rip event, can be avoided. Perhaps one of the most un- 
expected findings of the present research, was to show 
that the unusual behavior of the energy density asso- 
ciated with phantom-like matter (mimicked here by an 
appropriate NLED Lagrangean), in conjunction with the 
Randall-Sundrum brane effects, results in a model where 
the laws of gravity are simultaneously modified in the 
UV and in the IR limits (UV/IR mixing of scales). The 
above mixing is distinctive only of models that modify 
the matter part of thee Friedmann equation. 
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Avanzado de Cosmologia (lAC) collaboration. R. G.- 
S. acknowledges partial support from COFAA-IPN and 
EDI-IPN grants, and sip-ipn 20100610. T. G. acknowl- 
edges also the MES of Cuba for partial support of the 
research. 



VII. APPENDIX: DYNAMICAL SYSTEMS 

Here we include brief tips of how to apply the dynam- 
ical systems tools to situations of cosmological interest. 
In order to apply these tools one has to follow the steps 
enumerated below: 



tial equations (ODE), say:'' 

2. With the help of the chosen phase space variables, 
to build an autonomous system of ODE out of the 
original system of cosmological equations (t is the 
time-ordering variable, not necessarily the cosmic 
time): 



Notice that the RHS of these equations do not de- 
pend explicitly on r (that is the reason why the 
system is called autonomous). 

3. To identify the phase space spanned by the chosen 
variables {xi,X2, ■■■Xn), that is relevant to the cos- 
mological model under study. This amounts, basi- 
cally, to define the range of the phase space vari- 
ables that is appropriate to the problem at hand: 

5* — {{xi, X2, ■■■Xn) : bounds on the Xi-s} . 

4. Finding the equilibrium points of the autonomous 
system of ODE, amounts to solve the following sys- 
tem of algebraic equations on {xi,X2, ■■■Xn)- 

fi{xi,X2, ...Xn) = . 

5. Next one linearly expands the equations of the au- 
tonomous system of ODE in the neighborhood of 
the equilibrium points pk = Pk{xi,X2, ■■■Xn), k = 
l,2,...m:^ I. e., one replaces Xi ^ Xi + e^, where 
Ci are the small (linear) perturbations around the 
equilibrium points. Hence the system of ODE be- 
comes a system of linear equations to determine the 
evolution of the e^-s: 

§^/-i:(||),.-(»i.. 

otherwise, since fi — fi{p) — 0, then 
j 

where we have introduced the linearization or Ja- 
cobian matrix [M/] = dfi/dxj. 



1. To identify the phase space variables that allow 
writing the system of cosmological equations in the 
form of an autonomus system of ordinary differen- 



* There can be several different possible choices, however, not all 
of them allow for the minimum possible dimensionality of the 
phase space. 

^ In general the number of equilibrium points is different from the 
dimension of the phase space: m n. 
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6. The next step is to solve the secular equation to de- 
termine the eigenvalues Xi of the linearization ma- 
trix at the given equilibrium point p: 

det|M(p)^'-AC//| = 0, 

where [U^] is the unit matrix. 

7. Once the eigenvalues of the linearization around a 
given equilibrium point p have been computed, the 
evolution of the perturbations is given by 

n 
j 

where the amplitudes (cq)^ are constants of inte- 
gration. 

If all of the eigenvalues have negative real parts, the per- 
turbations decay with r, i. e., the equilibrium point is 



stable against linear perturabtions. The corresponding 
equilibrium point is said to be a future attractor. If at 
least one of the eigenvalues has positive real part, the 
perturbations grow with r so that these are not stable in 
the direction spanned by the given eigenvalue. Hence the 
point is said to be a saddle. The perturbations around 
a given equilibrium point are unstable, in other words 
the point is a past attractor - also acknowledged as a 
source point in the phase space -, if all of the eigenval- 
ues have positive real parts. Points whose linearization 
is characterized by complex eigenvalues are said to be 
spiral equilibrium points, and are commonly associated 
with oscillatory behavior of the corresponding solution. 
If at least one of the eigenvalues has a vanishing real 
part, the equilibrium point is said to be non- hyperbolic. 
In the latter case, in general, and unless some of the 
non-vanishing real parts of the eigenvalues are of oppo- 
site sign, one can not give conclusive arguments on the 
stability of the equilibrium point. Other techniques have 
to be applied. 
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